LETTER TO THE EDITOR 



The stress— energy operator 

Adam D Heifer 

Department of Mathematics, Mathematical Sciences Building, 
University of Missouri, Columbia, Missouri 65211, U.S.A. 



Abstract. We compute the stress-energy operator for a scalar linear 
quantum field in curved space-time, modulo c-numbers. For the 
associated Hamiltonian operators, even those generating evolution along 
timelike vector fields, we find that in general on locally Fock-like 
('Hadamard') representations: (a) The Hamiltonians cannot be self- 
adjoint operators; (b) The automorphisms of the field algebra generated 
by the evolution cannot be unitarily implemented; (c) The expectation 
values of the Hamiltonians are well-defined on a dense family of 
states; but (d) These expectation values are unbounded below, even 
for evolution along future-directed timelike vector fields and even on 
Hadamard states. These are all local, ultraviolet, effects. 
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In a relativistic quantum field theory, two sorts of evolution are important: that of the 
state vectors; and that of the field operators. The state vectors are required to evolve 
by unitary motions of the Hilbert space (except when reduction occurs) in order that 
probability be conserved. On the other hand, the evolution of the field operators is 
governed by relativistic field equations on space-time, and need not arise from unitary 
motions. 

To spell out what this means, we recall that a relativistic field theory is constructed 
by first forming an algebra A of fields, encoding the canonical commutation relations. 
This algebra is abstract in that it is not (yet) represented as a collection of operators 
on a physical Hilbert space; the choice of such a representation is the second step in 
the construction. Now evolution in space-time is represented by flow along a vector 
field, and this flow induces a family of automorphisms of A. But these automorphisms 
need not arise from unitary motions on a given representation of A. We show here 
that this failure occurs in general when (as is usual) the representation is taken to be 
'Hadamard' — that is, to have the same leading short-distance behavior as the Fock 
representation in Minkowski space. 

The phenomenon of non-unitarily implementable evolution has been known for 
some time, but has generally been thought possible only when unaccounted-for 
massless radiation might escape to infinity (the 'infrared catastrophe'; see Bjorken 
and Drell 1965, Jauch and Rohrlich 1975, Ashtekar and Narain 1981). Here it 
occurs as a local (ultraviolet) effect. Since the evolution is supposed to be generated 
by Hamiltonian operators, it is perhaps not surprising that the lack of unitary 
implementability is tied to a lack of existence of fully satisfactory Hamiltonians. These 
cannot exist as operators on reasonable domains in the physical Hilbert space. Their 
expectation values are well-defined on a dense family of states, but generically these 
expectation values are unbounded below. 

The present results should be contrasted with those of Brunetti et al (1995). 
These authors proved that the stress-energy operator exists as an operator-valued 
distribution on Hadamard representations. This means that / Ta^f"^ dvol exists as a 
self-adjoint operator when is taken to be a compactly supported test-function, 
smooth on space-time. However, the Hamiltonian operators are given by integrating 
Tah over Cauchy surfaces (equivalently, by taking f°-^ to be a kind of distribution), 
and so the work of these authors is not in confiict with the present result. Comments 
about the physical significance of the two results will be found at the end of this letter. 

The arguments of Brunetti et al are general but not explicitly computational, 
and they establish the existence of the stress-energy defined by the point -splitting 
procedure (cf. Wald 1994). It will turn out that our stress-energy's expectations 
verify Wald's key axioms (his (1) and (2)), and then Wald's work shows that our 
stress-energy is the same as that of Brunetti et al (modulo c-numbers). 

This letter is given to outlining the arguments in a simple case, that of an 
uncharged spinless field of mass m in a curved space-time {M^Qab)- Similar effects 
occur for fields responding to scalar and electromagnetic potentials in Minkowski 
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space; these will be treated elsewhere, along with many mathematical details. We use 
the by-now standard symplectic quantization, which is treated well in Wald's (1994) 
book; see also Ashtekar and Magnon-Ashtekar (1980), Birrell and Davies (1982), 
DeWitt (1983), Fulling (1989). See Rccd and Simon (1972-5) for general operator 
theory and Treves (1980) for pseudodiffcrcntial operators. Our conventions for space- 
time quantities are those of Penrose and Rindler (1984-6). We set c = ^ = 1. 
The field equation is 



We assume the space-time is oriented, time-oriented and globally hyperbolic. We 
also take the Cauchy surfaces to be compact; this is done only for simplicity: it avoids 
some operator-theoretic technicalities and rules out all infrared effects. We let F be 
the space of smooth classical solutions to (1); it is equipped with a symplectic form 
a;((/), •0) = - 0*dV'). See Wald (1994) for a precise definition of A. 

The representations of interest are mathematically analogous to the Fock 
representation. Let J : F — > F be a map preserving uj, satisfying = — 1 and 



is positive-definite. (In Minkowski space, one chooses J0 = i(0+ — 0_), where 0± 
are the positive- and negative-frequency parts of 0.) Then J makes F (completed 
with respect to the norm (2)) a complex Hilbert space F(J), which is the analog of 
the one-particle Hilbert space, f and the full Hilbert space may be constructed from 
this as usual. In fact, it is a little more convenient to start with F(— J); then the 
full Hilbert space 7Y(J) may be identified with a space of holomorphic (rather than 
antiholomorphic) functions on F(J), the elements of the symmetric tensor powers of 
F(— J) forming the coefficients of the power-series expansion. On such a function $(0), 
the action of the field operators is determined by oj{a, 0)$(0) = a;(Q;_, 0)$ + IVq.^^, 
where a± are the J-linear and J-antilinear parts of the classical solution a. It should 
be noted that H{J) contains only the holomorphic functions whose norm (in the 
appropriate sense) is finite. 

The choice of J is important. In general, different choices will not lead to unitarily 
equivalent representations; one has unitary equivalence iff their difference is Hilbert- 
Schmidt, that is, iff tr(Ji — J2)^ < oo. Kay and Wald (1991) have shown that there 
is a choice of J such that the two-point functions {(j){x)(l){y)) of the J-' vacuum' have 
the same asymptotic form as a; — > y (in geodesic normal coordinates) as for the 
Fock representation, and that (in the case of compact Cauchy surfaces) this choice 
is unique up to unitary equivalence. Such a representation is called Hadamard. For 



(V„V'^ + m2)(/. = 0. 



(1) 



such that 



a;((/), J0) 



(2) 



f This analogy is in general only mathematical. The identification of which mathematical structures 
correspond to physical one-particle states is involved and to some degree ill-defined. See Wald (1994) 
and references therein. 
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the remainder of this paper, we assume the physical Hilbert space is a Hadamard 
representation. 

We now take up the problem of evolution. The algebra A is generated by the 
Cauchy data for (p on any Cauchy surface. Flowing along a vector field in M induces 
a canonical transformation of Cauchy data, which generates an automorphism of A. 
Our first result is that in general this automorphism need not be induced by a unitary 
map on 7i(J). In other words, the automorphism may be regarded as determining a 
deformation of the representation. 

Consider a space-time with M = {(t, x) |teM, xgE}, where t is a global time 
function and for t <t_ (resp., t >t+ > t_) the metric is static of the form 

dt^ - h^h^x'^dx'' , (3) 

where /i^^ is a t-independent positive-definite three-metric on E. The vector field ^ 
is taken to be d/dt, and the Cauchy surfaces to be surfaces of constant t. We shall 
write Ts = {(</>, tt) G C°°(S) x C°°(E,A^)} for the space of Cauchy data (using the 
'momentum' tt = *d(p pulled back to E in place of the normal derivative). Then we 
have isomorphisms : Fs ^ F embedding these as classical solutions with data given 
at time t. 

Suppose that evolution were unitarily implement able. Let J be a Hadamard 
complex structure on F. An allowable choice of J in the regime t < t- (resp., t > 
in block form with respect to the decomposition above, is it±J±it2, where 

(m2 - A±)-'/'(dvol(/i±))- 

-(dvol(/i±))(m2 - A±)^/^ 



J± = 



(4) 



and A± is the Laplacian associated with /i^^^. Evolution from t_ to t+ is given 
by the automorphism g = itj^it} of F. Our hypothesis is that this is unitarily 
implementable, and this implies that the complex structure git_J-it}g^^ is Hadamard. 
But git_J-if}g''^ = it+J-it^- Therefore the structures 

are all Hadamard, and by equivalence of Hadamard structures, the difference between 
any two must be Hilbert-Schmidt. Taking the difference of the last two, and 
conjugating by it^ (which preserves the Hilbert-Schmidt class) implies that J_ — J+ 
must be Hilbert-Schmidt. However, it is easy to see that this fails in general. (For 
example, if h^j^ = Ch^j^ for a constant C, then tr( J_ — J+)^ is explicitly calculable in 
terms of the eigenvalues of the spatial Laplacian.) This is a contradiction, and so the 
evolution cannot be unitarily implemented. 

The foregoing is a rigorous proof that a certain class of space-times, vector fields 
and Cauchy surfaces has non-unitarily implementable evolution. But one expects 
the result to hold generically, since tr( J — gJg~^Y can be expressed as an integral of 
products of Green's functions and kernels for the two-point function, so the fact that 
it diverges should be stable. 
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Although the evolution is not unitarily implement able, one can still search for a 
generator, that is, an operator H{^,'E) such that 



[#(^,S),0] = -irV„0 onS 



(6) 



(and the usual formula for [H{^, E), Vn^], with Vn the normal derivative). Such an 
operator does exist, although in a limited sense. We adopt a standard index notation 
for the complex vector space r(J); we write 0" for an arbitrary element of this space, 
and ijaa for the inner product. Then the classical Hamiltonian operator E) on F 
may be decomposed into a J-linear part iB°'i3 and an antilinear part C'^p. We set 
Aaf3 = {-i/4:)rfaaC"f}. Then formally 



xi3 d d 



(7) 



generates evolution. Note that H{^, E) is equal to the classical Hamiltonian function 
with (f) replacing (f) and normal-ordered with respect to the senses of creation and 
annihilation defined by J. Therefore it arises from a stress-energy operator Tab which 
is the J-normal ordered transcription of the classical stress-energy. It can be shown 
from this that the expectation values ($|Ta;,|$) verify Wald's (1994) axioms. One sees 
that (^, E) is determined modulo c-numbers by the requirement that it generate 
evolution. 

For stronger results, we need a usable formula for J. A straightforward lengthy 
calculation shows that in normal coordinates (for the induced metric on E), the symbol 
of J is given by 



sym J (C) = 



syma |C| ^ 
-|Cl -sym a 



+ 



o{\cn 

o(ICI-^ioglCI) 



o(|C|-=^iog|CI) 
o(IC|-^) 



where 2symQ! = |C|~''^7''a" — CaCb and TVab is the second fundamental form of E. 

(The operator J is obtained to corresponding accuracy by putting — iVa for (a- See 
Treves 1980.) With this in hand, it is also clearest to invoke the isomorphism of F( J) 
with the +i eigenspace of J in F C to identify with the complex wave-function 
0CX5 defined as the configuration value of (1/2) (1 — iJ)(j).] Then with s = y/rn^ — A, 
in the case ^" = where is the unit normal to E, we find 



,{s + l.o.t.)0cxdE 



Aapr<P'' = (i/4) J {aPs" + l.o.t.)(/>e. ■ </>cxdE 

-B0CX = (/5S + l.O.t.)0cx 



(9) 



t It is possible for the projection onto the configuration values to have a finite-dimensional kernel. 
However, this does not happen generically, and in any event it is not hard to modify what follows to 
apply to this case. 



6 



Letter to the Editor 



where 'l.o.t.' means 'lower-order terms' in the sense of pseudo-differential operators 
(Treves 1980). 

Remarks: (a) A computation like that of (9), but keeping all terms shows that 
B°'p is formally self-adjoint, and that for everywhere timelike, the operator B°'p is a 
real elliptic pscudodifferential operator and has a canonical self-adjoint extension, (b) 
From (9), we see Aap is a bounded symmetric form, but that generically AapA = oo, 
so -ff (^, S) cannot be an operator from any domain containing the constant functions 
to 7i{J). (c) However H{^,T,) does exist as a quadratic form on, for example, the 
space of polynomials each of whose terms is a tensor of finite rank over the dual of the 
domain of B^a- Such a domain is dense in H{J). Then H{^, E) can be given meaning 
as an operator from this domain to a certain space of polynomials, not elements in 
H{J). In this sense (6) holds. 

Our next result, perhaps the most disturbing, is this: that even for ^" timelike 
future-pointing, generically the form H(^, E) is unbounded below. Take to be the 
unit normal to S, for simplicity. Then using (9) one can show that states of the form 
$(0) = 1 + bai3(l)"(l)'^ with 

6a/30"/ = -kA^pr{B''(l>f , (10) 

with k a positive constant, for —3 < z/ < —2, have finite norm but E) |$) = 

— oo. Now bai34>'^4>'^ is given by integrating the Cauchy data for against a certain 
distributional kernel. We may approximate this distribution by a sequence of smooth 
kernels, giving a sequence of states ^j{4>)- Each of these is a Hadamard state, that 
is, has two-point functions with singular parts given by the classical Hadamard 
expansion (Wald 1994). The states are uniformly bounded in norm and have 
($,|^(e,E)|$,)^-oo. 

Finally, for E a surface of constant time in Minkowski space, if is not a Killing 
vector, similar calculations show that still the Hamiltonian fails to exist as a self- 
adjoint operator, and still is unbounded below as a form. (Compare Epstein et al 
1965.) Thus what distinguishes special from general relativity is the ability to find E 
and ^" such that the Hamiltonians have satisfactory properties. This ability depends 
on the existence of hypersurfaces with symmetries. (The leading pathogenic terms 
cancel iff there is a solution to Killing's equation projected to E.) 

We now discuss the consequences of these results. The central question they raise 
is. How is energy exchanged between the quantum field and gravitational (or other) 
fields? 

Insofar as it is possible to accept the external-field approximation as a model of 
physics, and deferring the identification of a self-adjoint energy operator, the results 
uncovered here are simply curiosities: they do not affect the dynamics of the theory. 
However, this neglects interactions; the real world would be described by the model 
subject to perturbations. If (J TabC^dE**) were any sort of measure of the energy of a 
state, one would expect the perturbations to lead to the emission of arbitrarily large 
energies from the quantum field, an effect which is not observed. (This concern is 
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parallel to that which led to the second quantization of the electron field, of course.) 
Either one must explain why perturbations do not have this effect, or one must 
abandon the interpretation of the expectation value as a measure of energy. 

The non-existence of the Hamiltonians as self-adjoint operators means that, in 
arbitrarily weak gravitational fields, we do not have a satisfactory definition within 
conventional quantum theory of the energy as an observable. This leads to two 
questions. Is there a satisfactory definition, perhaps of another nature? And, If 
J Tab^"dS* is not an observable, what quantity is measured in experiments? 

Some progress on these questions ought to be possible by making explicit enough 
models of energy measurement. Recall that the problems uncovered here occur in 
Minkowski space, if the vector field is not Killing. But a real experiment, which 
occupies a finite space-time volume, cannot measure / Tab^"'dT,^ for a Killing field 
(since then is not compactly supported). What we measure in the laboratory is 
some self-adjoint operator which approximates J Tab^"'dT,^. A model of measurement 
should identify these approximants, and lead to an understanding of what bounds, if 
any, there are on observable energies. (Compare Ford and Roman 1995.) 

A more speculative possibility, which will be considered elsewhere, is that 
measurement of energy effects a change of representation, the new representation being 
determined so that the Hamiltonian exists as a semibounded self-adjoint operator. 

All that has been used so far is the external-field approximation. If we assume as 
well the semiclassical approximation. 



we have all the previous problems, but additionally the negative-energy states give rise 
to arbitrarily strong repulsive gravitational fields. A possible remedy is to use a space- 
time averaged stress-energy operator. Indeed, probably many workers on semiclassical 
gravity have regarded (11) as a short-hand, the right-hand side being more properly 
replaced by a space-time average. Recent intriguing 'quantum inequalities' suggest 
that this might remove the negatively infinite energies (Ford and Roman 1995, see 
also Flanagan and Wald 1996). In this context, the force of the present results is that, 
since the limit of the average over smaller space-time volumes does not exist, there is 
no ideal value to which averages tend, and the details of the averaging must be spelled 
out. It seems worthwhile pointing out that, to have reasonable causality properties, 
one would expect Tab{x) to be replaced by an average over the causal past J~{x); this 
would introduce a manifest time-asymmetry into the theory. 

The most striking applications of quantum field theory in curved space-time — the 
prediction that black holes should lose mass by Hawking radiation, and the inflationary 
cosmologies — depend on the semiclassical approximation. We must regard these as 
at least thrown into question by the present results. 




(11) 
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